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A homogeneous Bose gas is investigated at finite temperature using renormalization 
group techniques. A non-perturbative flow equation for the effective potential is derived 
using sharp and smooth cutoff functions. Numerical solutions of these equations show 
that the system undergoes a second order phase transition in accordance with universality 
arguments. We obtain the critical exponent v = 0.73 to leading order in the derivative 
expansion. 
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I. INTRODUCTION 



The remarkable achievement of Bose-Einstein condensation (BEC) of alkali atoms in magnetic 
traps [1-3] has created an enormous interest in the properties of dilute Bose gases. A very recent review 
on the trapped Bose gases can be found in Ref. M and an overview in Ref. S. However, it is useful 
to gain insight into the simpler problem of an interacting homogeneous Bose gas by applying modern 
methods from thermal field theory before attacking the full problem of atoms trapped in an external 
potential. 

The homogeneous Bose gas at zero temperature was intensively studied in the 1950s (^J^]. The 
properties of this system can be calculated in an expansion of powers of \J pa 3 , where p is the density of 
the gas and a is the S-wave scattering length. At zero temperature this expansion is equivalent to the 
loop expansion. The leading quantum correction to the ground state energy was calculated in 1957 by 
Lee and Yang The complete two-loop result was recently obtained by Braaten and Nieto ||. 

The system has also been investigated at finite temperature. Since this model is in the same univer- 
sality class as the three-dimensional xy-model, one expects the phase transition to be second order |J. 
However, both Bogoliubov theory Q and two-body t-matrix theory [ flO| predict a first order phase tran- 
sition. They fail because they do not take into account many-body effects of the medium |ll[ . In order to 
resolve this problem, Bijlsma and Stoof used a many-body t-matrix approximation fll|| . In this approx- 
imation the propagators are self-consistently determined in the self-energy diagrams (in contrast with 
the Bogoliubov theory, where free propagators are used in the self-energy graphs). This approach yields 
a second order phase transition, but predicts the same critical temperature as that of a non-interacting 
gas. 

Haugset, Haugerud and Ravndal jt2) have recently studied the phase transition of this system. By 
self-consistently solving a gap equation for the effective chemical potential, they are effectively summing 
up daisy and superdaisy diagrams. The inclusion of these diagrams is essential in order to satisfy the 
Goldstone theorem at finite temperature. Within this approximation the phase transition is second order, 
but there is no correction to T c compared with the non-interacting gas. 



One very powerful method of quantum field theory is the Wilson renormalization group (RG) [ |13|Jl4| . 
Renormalization group techniques have been applied to a homogeneous Bose gas by several authors [15- 
18]. Bijlsma and Stoof have made an extensive quantitative study of this system and in particular 
calculated non-universal quantities such as the critical temperature and the superfluid fraction below 
T c The non-universal quantities depend on the details of the interactions between the atoms. Using 
renormalization group methods, they demonstrate that the phase transition is indeed second order, and 
that the critical temperature increases by approximately 10% for typical alkali atoms compared with that 
of a free Bose gas. A review summarizing the current understanding of homogenous Bose gases can be 
found in |o|. 

The critical exponents for the phase transition from a superfluid to a normal fluid observed in liquid 
4 He have been measured very accurately 0. On the theoretical side, the most precise calculations to 
date involve the e-expansion. The agreement between the five loop calculations up to e 5 and experiment 
is excellent, but one should bear in mind that the series is actually asymptotic. The e-expansion works 
extremely well for scalar theories, but not for gauge theories 0], and so it is important to have alternative 
methods to compute the critical exponents. The momentum shell renormalization group provides one 
such alternative, and the literature on the calculational techniques for obtaining the critical exponents is 
now vast [22-26] . 
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In the present work we reconsider the nonrelativistic homogeneous spin zero-Bose gas at finite 
temperature using RG techniques including higher order operators than was done in Ref JTs(] . We focus 
in particular on the critical behavior and the calculation of critical exponents. The paper is organized 
as follows. In section II we briefly discuss the symmetries and interactions underlying the effective 
Lagrangian of the Bose gas. In section III the renormalization group is discussed, and we derive the 
flow equation for the one-loop effective potential. In section IV the flow equation for the one-loop RG- 
improved effective potential is found. In section V we calculate fixed points and critical exponents, using 
different cutoff functions and we address the question of scheme dependence of the results for the critical 
exponents j25j. In section VI we summarize and conclude. In the Appendix we prove that the one-loop 
renormalization group flow equation is exact to leading order in the derivative expansion. Throughout 
the paper we use the imaginary time formalism. 



II. EFFECTIVE LAGRANGIAN 



In this section we briefly summarize the symmetries and interactions underlying the effective La- 
grangian of a dilute Bose gas. A more detailed description can be found in [5|p7|]. 

The starting point of the description of a homogeneous Bose gas is an effective quantum field theory 
valid for low momenta p8[ | . As long as the momenta p of the atoms are small compared to their inverse 
size, the interactions are effectively local and we can describe them by a local quantum field theory. Since 
the momenta are assumed to be so small that the atoms are nonrelativistic, the Lagrangian is invariant 
under Galilean transformations. There is also an 0(2)-symmetry, and for simplicity we take the atoms to 
have either have zero spin or to be maximally polarized. We can then describe them by a single complex 
field: 

^ = -^=[^+^2]- (1) 



The interactions of two atoms can be described in terms of a two-body potential V(yi~Y2). The potential 
is repulsive at very short distances and there is an attractive well. Finally, there is a long-range Van der 
Waals tail that behaves like 1/R 6 . For the example of 87 Rb the minimum of the well is around Rq = 5ao, 
where ao is the Bohr radius. The depth is approximately 0.5 eV and there are around 100 molecular 
bound states in the well. The S-wave scattering lengths a of typical alkali atoms are much larger than 
Ro (e.g. a 110a for 87 Rb) because the natural scale for the scattering length is set by the Van der 
Waals interaction ||. 

The Euclidean effective Lagrangian then reads 

C E = ^d T tp + • W> - A»VV + gii'Hf + ■■■ ■ (2) 



Here, /x is the chemical potential. We have set h = 1 and 2m = 1. The interaction g(ip^ip) 2 represents 
2^2 scattering and the coupling constant g is proportional to the S—w&ve scattering length a: 

g = 4™. (3) 



The ellipses indicate all symmetry preserving operators that are higher order in the number of fields and 
derivatives. 
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In the following we consider the dilute gas pa 3 <§; 1, which implies that we only need to retain the 
quartic interaction in the bare Lagrangian Eq. (||) [jl9| . 

The action can be written as 

S[ipi,M = £ ' drj d d x\^e ij ^ i d T i 3j + ±Wi ■ + f (Vi^i) 2 j, (4) 

where d is the number of spatial dimensions and repeated indices are summed over. 

In a field theoretic language, BEC is described as spontaneous symmetry breaking of the 0(2)- 
symmetry and the complex field tp acquires a nonzero vacuum expectation value v. Due to the 0(2)- 
symmetry, the field v can be chosen to be real and so we shift the fields according to 

ipi -> v + ipi, ip2^ip2- (5) 

Inserting (||) into (Q) and dividing the action into a free piece S{ TCC [ipi, ^2] and an interacting part 
•Sint [ipi , 1P2] we obtain 

5 free [u,^i,^2] = £dr J d d x i-\»v 2 + ^itijdripj + [-V 2 + V + V"v 2 ] ^ 

+ ^ 2 [-v 2 + y']^ 2 | (6) 

Sint[v, M2] = J 1 dr J d d x\ I [v 4 + Av^\ + 4^ 2 3 + + + ^1] 1 . (7) 



Here V(v) is the classical potential 



V(v) = -^v 2 + 9 l v\ (8) 



We will use primes to denote differentiation with respect to v 2 /2, so that V' = —fi + gv 2 and V" = 2gv 2 . 
The free propagator corresponding to Eq. (^) is a 2 x 2 matrix and in momentum space it reads 

where 

e P =P 2 
LU n = 2~KnT 



uJ 



V 



^{e P + V' + V"v 2 ][e p + V']. (10) 



In the broken phase the quadratic part of the action, Eq. (|6|) , describes the propagation of Bogoliubov 
modes with the dispersion relation 
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lu p = p^ep + 2fi. 



(11) 



The dispersion is linear in the long wavelength limit, corresponding to the massless Goldstone mode 
(phonons). This reflects the spontaneous symmetry breakdown of the 0(2) symmetry. In the short 
wavelength limit the dispersion relation is quadratic and that of a free nonrelativistic particle. 



III. THE ONE-LOOP EFFECTIVE POTENTIAL AT FINITE TEMPERATURE 



We are now ready to calculate quantum corrections to the classical potential Eq. (10). In this section 
we compute the one-loop effective potential which we will "RG improve" in the next section. This method 
of deriving RG flow equations is conceptually and technically simpler than the direct application of exact 
or momentum-shell RG techniques which is demonstrated in the Appendix. 

The one-loop effective potential reads 

Uf}(v) = V(v) + TrkiA- 1 (Lj n ,p) 

1.. 2 , 9 4 , l^V^ f dd P , r 2 2i 



2 -^ + |»« + jT^y^.„K +4 (i2) 



The sum is over the Matsubara frequencies, which take on the values uj n — 2imT ', and the integration is 
over d-dimensional momentum space. 

We proceed by dividing the modes in the path integral into slow and fast modes separated by an 
infrared cutoff k. This is done by introducing a cutoff function Rk{p) which we keep general for the 
moment. By adding a term to the action Eq. (Q): 

Sp, k [il>i,ih] = S[ih,ih] + J dr J d 3 z±J4(v /= V2)V^ • (13) 

the modified propagator reads 

a / s 1 / e p (R k (p) + 1) + V uj n \ 

AfcK ' P) = 1%+* { e p (R k ( P ) + 1) + V + V»v* ) < (U) 

and the modified dispersion relation is 

u P ,k - \J[t P {Rk{p) + 1) + V' + V"v^\ [e p (R k (p) + 1) + V']. (15) 



By a judicious choice of Rk(p), we can suppress the low momentum modes in the path integral and leave 
the high momentum modes essentially unchanged. In section |v| we return to the actual choice of cutoff 
functions. It is useful to introduce a blocking function fk(p) which is defined through 

*W - i^m. da) 
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The blocking function satisfies 



lim fk(p) = 0, 



lim f k {p) = 1. 

p—*OG 



(17) 



This implies that the function Rk(p) satisfies 

lim Rk(p) = co, 

p^O 



lim R k {p) = 0. 

p—>OG 



(18) 



These properties ensure that the low-momentum modes are suppressed by making them very heavy and 
the high-momentum modes are left essentially unchanged. 



We return to the one-loop effective potential, which now becomes 



(19) 



Here, the subscript fc indicates that the effective potential depends on the infrared cutoff. Upon summa- 
tion over the Matsubara frequencies, we obtain 



UpAv) = V(v) + 



1 f d d p 



(2it) d 



j p , fc + 2Tm[l- p-^-x] 



(20) 



The first term in the brackets is the T = piece and represents the zero-point fluctuations. The second 
term includes thermal effects. Differentiation with respect to the infrared cutoff k yields: 



k f d d p (dR k 



(27r) d V dk 



1 



1 



2u) p ,k uJpAe 13 ^-" - 1) 



[2e p , fc + 2V' + V"v 2 ] 



(21) 



Eq. ( pi| ) is the differential equation for the one-loop effective potential. It is obtained by integrating 
out each mode independently, where the feedback from the fast modes to the slow modes is completely 
ignored. Since all modes are integrated out independently, this is sometimes called the independent mode 
approximation p9j. 



Equation (20) provides an inadequate description of the system at finite temperature in several 
ways. Since the minimum of the one-loop effective potential at finite temperature is shifted away from 
the classical minimum, the Goldstone theorem is not satisfied. This theorem is known to be satisfied for 
temperatures below T c to all orders in perturbation theory ]3C| , and any reasonable approximation must 
incorporate that fact. 

Secondly, it is clear from Eqs. ( |l5| ) and ( pp| ) that the one-loop effective potential has an imaginary 
part for all temperatures and for sufficiently small values of the field v, when the bare chemical potential 
is positive. However, we know that a thermodynamically stable state for T > T c corresponds to v = 
and so the effective potential is purely real for sufficiently high temperatures. More generally, ordinary 
perturbation theory breaks down at high temperature due to infrared divergences and this has been 
known since the work on summation of ring diagrams in nonrelativistic QED in 1957 by Gell-Mann and 
Bruckner In the next section we derive an RG equation, whose solution has none of the above 

shortcomings. 
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IV. RENORMALIZATION GROUP IMPROVEMENT 



In the previous section we derived the one-loop effective potential at finite temperature and discussed 
the fact that it is not capable of reliably describing the system at finite temperature. The lack of feedback 
from the fast modes to the slow modes as we lower the infrared cutoff k leads to a poor tracking of the 
effective degrees of freedom causing the problems mentioned above. The situation is remedied by applying 
the renormalization group, which effectively sums up important classes of Feynman diagrams |26|. In 
order to obtain the differential equation for the RG-improved effective potential, we do not integrate out 
all the modes between p = oo and p = k in one step. Instead, we divide the integration volume into small 
shells of thickness Ak, then lower the cutoff from k to k — Ak and repeat the one-loop calculation. This 
is equivalent to replacing V by Up t k on the right hand side of Eq. (pi]), making it self-consistent pi: 



,d __k f d d p (dR k 
k dk UfJ - k ~ 2 J (27T)* ^ dk 



1 



1 



2uj Ptk uj p . k (e^p- k - 1) 



[2e p 



2U kk 



(22) 



where 



J p,k 



i p {R k {p) + l) 



,(R k (p) + l) 



(23) 



and the primes in Eqs. ( |22| ) and ( p3| ) denote differentiation with respect to v 2 /2. The self-consistent 
Eq. ( p2[ ) is not a perturbative approximation, but is exact to leading order in the derivative expansion. 
This equation is derived in the Appendix without performing a loop expansion. 



Note that since 



2v- 



dU, 



0,k 



dU, 



dv 2 



dv 



(24) 



the dispersion relation at the minimum of the effective potential in the broken phase reduces to 

u p ,k=o = PyJe P + U p,k=o v2 - ( 25 ) 



Hence, the Goldstone theorem is automatically satisfied for temperatures below T c . 

This equation interpolates between the bare theory for k = oo and T — and the physical theory at 
temperature T, for k = 0, since we integrate out both quantum and thermal modes as we lower the cutoff. 
This implies that the boundary condition for the RG-equation is the bare potential, Up^ooW) — V(v). 

In Refs. J3^,|33| renormalization group ideas have been applied to A0 4 theory using the real time 
formalism. In the real time formalism one can separate the free propagator into a quantum and a thermal 
part |Q, and in [^2|,^3| the infrared cutoff is imposed only on the thermal part of the propagator. This 
implies that the theory interpolates between the physical theory at T — and the physical theory at 
T ^ 0. Hence, the boundary condition for the RG-equation in this approach is the physical effective 
potential at T = 0. However, if one imposes the infrared cutoff on the both quantum and thermal part 



of the propagator, one can derive Eq. (22), showing that identical results are obtained using the two 
formalisms. 
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We close this section by commenting on the choice of cutoff function. It is clear from Eq. (|2^) that 
the non-perturbative flow equation depends explicitly on the choice of Rk (p) ■ We know that the nonzero 
Matsubara modes are strongly suppressed at high temperature and can be integrated out perturbatively; 
the important point is to treat the zero mode correctly. For a thorough discussion of various finite 
temperature cutoff functions applied to relativistic A0 4 theory see Ref. [p5[ . 



V. RESULTS 



In this section we present our results for the numerical solution of the renormalization group flow 
equation and the calculations of the fixed point and critical exponents. We consider the cases of a sharp 
cutoff and a smooth cutoff separately. 



A. Sharp Cutoff 



The sharp cutoff function is defined by the blocking function fk(p) = Q{p — k), which is displayed in 
Fig [l] (solid line). It provides a sharp separation between fast and slow modes. Using the sharp cutoff 
the slow modes become completely suppressed in the path integral, while the fast modes are completely 
unaltered. The adva ntag e of using the sharp cutoff function compared to the smooth cutoff functions 
considered in section VE is that the integral over p can be done analytically, resulting in a differential 
RG-equation. In this case Eq. (B2h reduces to 



k dk U ^ k 



S d k d 



2Tln[f 



(26) 



Here, 



S d = 



(47r) d / 2 r(d/2) ' 



(27) 



Eq. (H^) is derived in the Appendix. 

We have solved Eq. (|26| ) numerically for d — 3 and the result for different values of T are shown 
m Fig | The curves clearly show that the phase transition is second order. For T < T c , the effective 
potential has a small imaginary part, and we have shown only the real part in Fig. |^. The imaginary 
part of the effective potential does, however, vanish for T > T c in contrast to the independent mode 
approximation in which it does not. The effective chemical potential Hf),k as well as the quartic coupling 
constant gp t k (defined as the discrete first and second derivatives of the effective potential with respect to 
v 2 /2) are displayed in Fig. [}]and both quantities vanish at the critical point. The corresponding operators 
are relevant and must therefore vanish at T c , and we see that the renormalization group approach correctly 
describes the behavior near criticality. Moreover, the sectic coupling goes to a non-zero constant 
at T c . The inclusion of wavefunction renormalization effects turns the marginal operator gjp k into an 
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irrelevant operator that diverges at the critical temperature pq] . The success of describing the phase 
transition using the renormalization group is due to its ability to properly track the relevant degrees of 
freedom. The dressing of the coupling constants as we integrate out the fast modes is taken care of by 
the renormalization group and this is exactly where the independent mode approximation fails. 

In order to investigate the critical behavior near fixed points, we write the flow equation in dimen- 
sionless form using 

P = (3k 2 
v = f3 1 ' 2 k d ~ 2 v 
% fe = 0k- d Uf,, k 
uj k = k~ 2 uj k . (28) 



This yields 



dk 



(29) 



The critical potential is obtained by neglecting the derivative with respect to k on the left hand side of 
Eq. (p9|). Expanding in powers of /3u)fc we get 



s, 



u p,k = — irfiuk - S d \n [/3, 



(30) 



Taking the limit /3 — s- and ignoring the term which is independent of v leads to 

Sd, 



u P,k- y 



In [l + U'] + ln [l + U'+ U"v 2 ] 



(31) 



This is exactly the same equation as obtained by Morris for a relativistic 0(2)-symmetric scalar theory 
in d dimensions to leading order in the derivative expansion | f36|| . Therefore, the results for the critical 
behavior at leading order in the derivative expansion will be the same as those obtained in the d— 
dimensional 0(2)-model at zero temperature. 

The above also demonstrates that the system behaves as a e?-dimensional one as the temperature 
becomes much higher than any other scale in the problem (dimensional crossover). This is the usual 
dimensional reduction of field theories at high temperatures, in which the nonzero Matsubara modes 
decouple and the system can be described in terms of an effective field theory for the n = mode in 
d dimensions p^] . The effects of the nonzero Matsubara modes are encoded in the coefficients of the 
three-dimensional effective theory. 

The RG-equation ( pl| ) satisfied by Up.k[v\ is highly nonlinear and a direct measurement of the 
critical exponents from the numerical solutions is very time-consuming. This becomes even worse as one 
goes to higher orders in the derivative expansion and so it is important to have an additional reliable 
approximation scheme for calculating critical exponents. In the following we perform a polynomial 
expansion |3q] of the effective potential, expand around v = 0, and truncate at ./Vth order: 



1 



(32) 
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The polynomial expansion turns the partial differential equation (26) into a set of coupled ordinary 
differential equations. In order to demonstrate the procedure we will show how the fixed points and 
critical exponents are calculated at the lowest nontrivial order of truncation (N — 2). We write the 
equations in dimensionless form using Eq. (pq) and 



9p,k = P^^Qp.k- 



(33) 



We then obtain the following set of equations: 
d _ 



2(l-^, fc ) 



[2n(Q k ) + 1] + Pn[Q k ) [n{Q k ) + 1] 



(34) 



Here, e = 4 — d and n(oj k ) is the Bose-Einstein distribution function written in terms of dimensionless 
variables 



n(wfe) 



1 



(35) 



A similar set of equations has been obtained in Ref. |19[ ] by considering the one-loop diagrams that 
contribute to the running of the different vertices. They use the operator formalism and normal ordering 
so the zero temperature part of the tadpole vanishes. 



The equations for the fixed points are 

d__ 

''dk P P< k ~~"' "dk 



k J^.H,k = 0, k—gp tk = 0. (36) 



Expanding in powers of (3{1 — fJ,g jk ) one obtains 



- -2 

o- 9p, k 1 „ 9p,k 5 

2 ^-^T^-r°> ^-^(t^f =0 - (37) 



If we introduce the variables r and s through the relations 



the RG-equations can be written as 



(38) 



— = -2 [1 + r] [r - S d s] , — = -s [e + 4r - 9S d s] . (39) 

We have the trivial Gaussian fixed point (r, s) = (0, 0) as well as the infinite temperature Gaussian fixed 
point (—1, 0). Finally, for e > there is the infrared Wilson-Fisher fixed point (e/5, e/ (5Sd)) ■ 
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Setting e = 1 and linearizing Eq. ( |39| ) around the fixed point, we find the eigenvalues (Ai,A2) = 
(—1.278, 1.878). The critical exponent v is given by the inverse of the largest eigenvalue; v = I/A2 = 0.532. 
This procedure can now be repeated including a larger number, N, of terms in the expansion Eq. j3^). 
The result for v is plotted in Fig. |] as a function of the number of terms, N, in the expansion. Our result 
agrees with that of Morris, who considered the relativistic 0(2)-model in d = 3 at zero temperature ]36| . 
The critical exponent v oscillates around the average value 0.73. The value of v never actually converges 
as N — - > 00, but continues to fluctuate. As Morris has pointed out in the .^-symmetric case, these 
oscillations are due to the presence of a pole in the complex v plane in the corresponding fixed point RG 
equation |^,^3| . Our results should be compared to experiment ( 4 He) and the e-expansion which both 
give a value of 0.67 Q. One expects that the critical exponent v converges towards 0.67 as one includes 
more terms in the derivative expansion. 



B. Smooth Cutoff 



In the previous section we considered the sharp cutoff function that divided the modes in the path 
integral sharply into slow and fast modes separated by the infrared cutoff k. However, there are alternative 
ways of doing this. In this section we consider a class of smooth cutoff functions -R™(p) defined through 

fi n iP) = (40) 

In the limit m — > 00 we recover the sharp cutoff function. A typical smooth blocking function is shown in 
Fig. (dashed line). We see that the suppression of the slow modes is complete for p = and gradually 
decreases as we approach the infrared cutoff. Similarly, the high momentum modes are left unchanged 
for p = 00 and there is an increasing suppression, albeit small, as one approaches k. Since we cannot 
carry out the integration over p analytically in Eq. (^2|), the RG flow equation is now more complicated. 
Taking the limit j3 — > and making a polynomial expansion as in the preceding subsection, we obtain 
the following set of dimcnsionlcss equations for N = 2: 



k iz9p, k = -eg + -^Sr [ Jl + ^H,k] ■ ( 41 ) 



dk 



where 



g 3 (s, m)s n ds 



l/m 



U ^ k) = Jo [s-,% k + 9 Hs,m)]^ ^ - {—) ^ 

In the case of a smooth cutoff function, we cannot calculate the fixed points and critical exponents 
analytically, but have to resort to numerical techniques. In Fig. |^ we have plotted the m-dependence 
of v for different truncations. Note in particular the strong dependence of m for N = 10. In Fig. ^. 
we have displayed the critical exponent v as a function of the number of terms N in the polynomial 
expansion using a smooth cutoff with m — 5 (solid line). For comparison we have also plotted the result 
in the case of a sharp cutoff (dashed line). The value of v continues to fluctuate, but the oscillations are 
significantly smaller for the smooth cutoff, and the convergence to its asymptotic range is much faster. 
Again, one expects the value of v to converge to the value 0.67 as more terms in the derivative expansion 
are included. 
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VI. SUMMARY AND DISCUSSIONS 



In the present paper we have applied renormalization group methods to the nonrelativistic homoge- 
neous Bose gas at finite temperature. We have explictly shown that the renormalization group improved 
effective potential does not suffer from the two major flaws of the one-loop effective potential; the Gold- 
stone theorem is automatically satisfied and the effective potential is purely real for temperatures above 
T c . The second order nature of the phase transition and the vanishing of the relevant couplings at the 
critical temperature have also been verified numerically. 

Truncating the RG equations at leading order in the derivative expansion, we have investigated 
the critical exponent v as a function of the number of terms N in the polynomial expansion of the 
effective potential and the smoothness of the cutoff function. In particular, we have demonstrated that 
the oscillations around the value v = 0.73 depends on the smoothness of the cutoff function, and that the 
oscillatory behavior can be improved by appropriately choosing the smoothness. The value m = 5 seems 
to be the optimal choice among the smooth cutoff functions investigated in the present paper. Whether 
the dependence on m is reduced as one goes to higher orders in the derivative expansion is not clear at 
this point. 

It is important to point out that it is not sufficient, as is conventional wisdom, to include only the 
relevant operators and perhaps marginal ones when calculating the d = 3 exponents. Instead, one has to 
make a careful study of the convergence of the exponents in question, as we have demonstrated. 

The present work can be extended in several ways. Expanding around the minimum of the RG 
improved effective potential instead of the origin is one posibility. This has been carried out in Ref. f23[| 
in the ^-symmetric case and the rate of convergence as function of N is larger. However, in the O(N)- 
symmetric case this expansion is complicated by the presence of infrared divergences due to the Goldstone 
modes |33| , and at present we do not know how to address that problem (sec also jl9| ) . 

The inclusion of wave function renormalization effects by going to second order in the derivative 
expansion will close the gap between the critical exponents of experiment and the e-expansion on one 
hand and the momentum shell renormalization group approach on the other. It is also of interest to 
investigate the influence of these effects on nonuniversal quantities such as the critical temperature and 
the superfluid fraction in the broken phase. One can also study finite size effects by not integrating down 
to k = 0, but to some k > where 1/k characterizes the length scale of the system under consideration. 
Of course, the real challenge is to describe the trapped Bose gas using renormalization group techniques. 



ACKNOWLEDGMENTS 



The authors would like to thank E. Braaten and S.-B. Liao for useful discussions. This work was 
supported in part by the U. S. Department of Energy, Division of High Energy Physics, under Grant 
DE-FG02-91-ER40690, by the National Science Foundation under Grants No. PHY-9511923 and PHY- 
9258270, and by a Faculty Development Grant from the Physics Department of The Ohio State University. 
J. O. A. was also supported in part by a NATO Science Fellowship from the Norwegian Research Council 
(project 124282/410). 



12 



APPENDIX: 



In this Appendix we give a proof that the renormalization group equation (|26|) is exact to leading 
order in the derivative expansion. 



The path integral representation of the partition function is 



(Al) 



where we have modified the action by adding a term to the action according to Eq. ([13]). The function 
G/3,k[j] is the generator of connected diagrams. Taking the derivative with respect to k 7 using Eq. (Al) 
and going to momentum space, we find that Gp tk [j\ satisfies the differential equation 



d d p dR k (p) j SG Pik [j]dG , k [j} 
(2ir) d dk \ 6ji(p) Sji(-p) 



Tr 



Sji(p)Sjj(-p) 



(A2) 



The symbol Tr indicates taking the trace over internal indices. 

The effective action T^^^] is defined through the Legendre transform: 

V{ = (ipi) 
= 6Gf, lk \j] 

Tp, k [v] = Gp tk \j] - J dr J d d x 2l v l - J dr J d d x \R k (^f^)V Vl • Vd, 



(A3) 
(A4) 



The last term in Eq. (A4) removes the additional term in Eq. ( |l3| ) from the effective action. The flow 
equation for T^^^] is 



9 r r l L tV f d P 



d d p dR k {p) 



dk 



e p Tr 



R k (p)e p Sij 



Svi(p)Svj(-p) 



(A5) 



To proceed we employ the derivative expansion of the effective action r^^fu]: 

Tf>, k [v] = ( dr [ d d xiu , k (v) + "-Z^ajVidrVj + \zfl(v)(V Vi f + 



(A6) 



The leading order in the derivative expansion is defined by setting the coefficients Zp\(v) and zj^ k (v) to 
unity and the coefficients of all higher derivative operators to zero. We denote the matrix in the brackets 
in Eq. (Al) by M and it reads 



M 



e P (R k (p) + l) + U'^ k + Ul k v 2 



e p (R k (p) + l) + U' 0ik 



(A7) 
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We have used the 0(2) symmetry to rotate v so that it points along the x-axis. Using the fact that 

Rkip) + 1 ' 



fk( P y 



the trace of M 1 can be written as 



TtM- 1 = 



f k {p)[2e p + 2U' p J k {p) + Ul k v 2 f k (p)] 



*v + U'^ k h{p)][e P + U' p J k {p) + Ul k v*ft(p)] + w*f k (p)S 



(A8) 



This yields 



d d p 



df k (p) 



{2-kY fkip) V dk 

2e P + 2Ui f3k f k (p) + U>>^f k (p) 

£p + U'^fkipWv + U'pjkip) + U'^Mp)} + ulflip) 



(A9) 



By doing the Matsubara sum, we obtain Eq. (p2|). If f k (p) is a function of only the ratio p/k then 

u df k (p) df k ( P ) 

K — = — P - 



Ok 



dp 



(A10) 



A change of variables t = f k (p) and using e p (t) = p 2 (t) yields 



dt 

T 



2p\t)+2U' 0ik t + Ul k vH 



lP 2 (t) + U^ k t]\p>{t) + U'^t + U« >k vH] + <4 



t 2 



(All) 



In the sharp cutoff limit, f k {p) — > 6(p—k), p{t) — > fc and performing the integral over £, Eq. (All) reduces 
to 



(A12) 



Summing over the Matsubara frequencies yields Eq. (& 
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FIGURE CAPTIONS: 

FIG. 1. The sharp blocking lunction (solid line) and a typical smooth blocking function (dashed line). 

FIG. 2. The real part of the RG-improved effective potential U@,k=o(v) for different values of the temperature. 
The phase transition is clearly second order. 

FIG. 3. The effective chemical potential Hp,k=o and the effective quartic coupling gp,k=o near the critical 
temperature. Both vanish at T c . 

FIG. 4. The critical exponent v as a function of number of terms, N, in the polynomial expansion. 

FIG. 5. The critical exponent v as a function of the smoothing parameter m for different values of the number 
of terms, N, in the polynomial expansion. 

FIG. 6. The critical exponent v as function of the number of terms , N, in the polynomial expansion using a 
sharp cutoff and a smooth cutoff with m = 5. 
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FIG. 1. The sharp blocking function (solid line) and a typical smooth blocking function (dashed line). 
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FIG. 3. The effective chemical potential fj,p t k=o and the effective quartic coupling gp,k=o near the critical 
temperature. Both vanish at T c . 
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. 4. The critical exponent v as a function of number of terms, N, in the polynomial expansion. 
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FIG. 5. The critical exponent v as a function of the smoothing parameter m for different values of the number 
of terms, N, in the polynomial expansion. 
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FIG. 6. The critical exponent v as function of the number of terms, N, in the polynomial expansion using a 
sharp cutoff and a smooth cutoff with m = 5. 
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